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ABSTRACT 

There is reason to suspect that about half of the baryons are in pressure-supported plasma in 
the halos of normal galaxies, drawn in by gravity along with about half of the dark matter. To be 
consistent with the observations this baryonic component, the corona, would have to be hotter 
than the kinetic temperature of the dark matter in the halo so as to produce acceptable central 
electron densities. We ascribe this hotter plasma temperature to the addition of entropy prior 
to and during assembly of the system, in an analogy to cluster formation. The plasma cooling 
time would be longer than the gravitational collapse time but, in the inner parts, shorter than 
the Hubble time, making the corona thermally unstable to the formation of a cloudy structure 
that may be in line with what is indicated by quasar absorption line systems. The corona of an 
isolated spiral galaxy would be a source of soft X-ray and recombination radiation, adding to the 
more commonly discussed effects of stars and supernovae. In this picture the mass in the corona 
is much larger than the mass in condensed baryons in a spiral galaxy. The corona thus would be 
a substantial reservoir of diffuse baryons that arc settling and adding to the mass in interstellar 
matter and stars, so that star formation in isolated spirals will continue well beyond the present 
epoch. 



Subject headings: cosmology; galaxies: halos 

1. Introduction 

Dynamics and weak lensing suggest that a sig- 
nificant fraction of the dark matter is in the viri- 
alized halos of i ^ i» galaxies. Gravity would 
have tended to draw a like fraction of the barons 
into the galaxies. Since the observed mass in 
condensed baryons — stars and the interstellar 
medium — is only about 8% of the total we have 
been led to consider the idea that roughly half of 
the baryons are in pressure-supported plasma in 
the outer parts of the galaxies (Fukugita, Hogan & 
Peebles 1998; Fukugita & Peebles 2004, FP). This 
is a modern variant of Spitzer's (1956) galactic 



corona, and the extent of the coronae we are con- 
sidering, about 300 kpc radius, is in line with the 
proposal by Bahcall & Spitzer (1969) that galac- 
tic coronae may be responsible for quasar absorp- 
tion lines. The massive corona picture is related 
to the evidence that X-ray detected groups con- 
tain cosmologically appreciable mass in plasma 
(Mulchaey 2000), and the search for a possibly 
significant plasma mass in the Local Group (Ras- 
mussen & Pedersen 2001 and references therein). 
Since groups are not much larger than the virial- 
ized parts of an individual galaxy the constraints 
on the plasma properties are similar. Here we are 
considering the possibility that there is a cosmo- 
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logically significant mass in the coronac of indi- 
vidual L ^ Lt, spiral galaxies, and our goal is 
to present an observationally acceptable model for 
the present-day structure of such a massive; corona. 

The starting idea for the model is that the bulk 
of the baryonic mass of a galaxy is in plasma gravi- 
tationally bound to the galaxy at near hydrostatic 
pressure equilibrium and at density low enough 
that the cooling time is longer than the Hubble 
time. At smaller radii the shorter cooling time 
would cause some of the plasma to settle toward 
the center of the galaxy, tending to relax to a near 
steady flow, in analogy to a cooling flow of intra- 
cluster plasma (e.g. Fabian 1988; Sarazin 1988). 
The first point to consider is that the plasma cool- 
ing time for most of the baryons in a virialized halo 
can be longer than the Hubble time, meaning the 
massive coronae can still be present if initial con- 
ditions allow the needed minimal clumping of the 
plasma. This is discussed in §2. 

Our model for the halo mass distribution to 
which the corona is gravitationally bound diff'ers 

from FP in two ways. First, we adjust the velocity 
dispersion to the value appropriate for an L ~ 
spiral galaxy. Second, we allow for the possibility 
of a less sharp truncation of the edge of the halo 
mass distribution, which we show can produce a 
more reasonable-looking galaxy-mass cross corre- 
lation function. We allow the possibility that heat- 
ing of the corona as it collapsed made it warmer 
than the temperature characteristic of the dark 
matter, in analogy to what is observed in clusters 
of galaxies (Allen et al. 2004). We assume that 
in the inner parts of the corona, where the cool- 
ing time is less than the Hubble time, the plasma 
typically has relaxed to a near steady state, not se- 
riously disturbed by recent merging events or by 
energy soiirces such as siipernovac, massive stars, 
or an active galactic nucleus. Since quasar absorp- 
tion lines show quite modest mass in atomic hy- 
drogen at radii > 30 kpc around galaxies (Chen et 
al. 2001), our model requires that the rate of con- 
version of plasma to condensed forms of baryons is 
significant only near the optically luminous parts 
of the galaxy. The model should also meet the re- 
quirement that the massive halo is consistent with 
soft X-ray observations of nearby galaxies and the 
integrated background. We construct two exam- 
ples of corona models under these conditions on 
the simple plan of spherically symmetric quasi- 



isothermal structures. Section 3 explains the mod- 
els and compares them to the observations. We 
discuss possible directions for further work in §4. 

There are alternative pictures for the present 
state of the baryons originally associated with the 
dark matter now in halos of galaxies. Shocks may 
have made most of this baryonic matter too hot 
to fall into massive halos. This is in line with 
the idea currently under discussion that there is a 
substantial baryon mass in warm-hot intergalactic 
matter at temperatures in the range 10^ — 10^ K 
(e.g. Dave et al. 2001), but it is difficult to see 
what could have produced a local shock power- 
ful enough to have prevented ongoing accretion of 
baryons by the Milky Way and Andromeda galaxy: 
the Virgo Cluster seems to be too far away, the 
galaxy velocity dispersion near the local sheet of 
galaxies seems to be too small, and if gravity pro- 
duced the local void it would not have been ac- 
companied by production of a shock. Other ideas 
assume that the baryons did fall into halos along 
with the dark matter. Possible, though we suspect 
unlikely, is the idea that most of these baryons col- 
lapsed to a condensed form that is not readily ob- 
servable, such as brown dwarfs. These condensed 
baryons would have to be distributed through the 
dark halo, because there is not room for them 
in the luminous parts of normal galaxies. This 
could be indicated by the extended diffuse glow 
of fight around M87 (Arp and Bertola 1969), but 
it requires a considerable departure from observed 
initial stellar mass functions; it may be related 
to the local MACHO population (Alcock et al. 
2000), though the indicated mass range, between 
0.15 and 0.9 Mq if a halo component, is a prob- 
lem. Silk (2003) develops the idea that supernova- 
driven winds may have blown most of the baryons 
out of the massive halos of protogalaxies. The en- 
ergy required is modest but, as Silk discusses, the 
efficiency of entrainment might be questioned. A 
low density wind would tend to drive convection 
through the coronae, leaving behind most of the 
baryons, while loc;ally deposited energy likely is 
radiated away rather than being adding to the en- 
tropy of the ambient material. Our conclusion is 
that the alternative, the massive corona picture, 
is worth exploring. 
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2. The Cooling Rate in the Outer Corona 

This discussion of the characteristic plasma 
cooling time for baryons in a virialized halo takes 
all parameters from FP. It is not significantly af- 
fected by the adjustments to the halo mass model 
discussed in the next section. 

We start with a measure of the plasma density 
at the nominal virial radius r2oo of a newly relaxed 
halo. The density of mass contained by r2oo is 200 
times the critical Einstein-de Sitter value, and the 
mass density p2oo at r2oo is one third of that if 
the density run within r2oo is p ^ r~^. To sim- 
plify this discussion we take it that the baryons 
and dark matter have the same spatial distribu- 
tion, so the mass fraction in baryons at r2oo is the 
ratio of the mass density parameters fife = 0.045 
and Qm = 0.28 in baryons and in all matter. (The 
baryon segregation discussed in §3.3 slightly re- 
duces the cooling time at r2oo-) With primeval 
helium abundance Yp = 0.248 the mean mass per 
particle is /i = 0.59 times the proton mass nip, 
and this measure of the plasma particle number 
density at r2oo is 



^200 



,limp 



1.0 X 10 



-4 cm-3, 



(1) 



with electron density ne,200 = 0.52n2oo- The one- 
dimensional mass- weighted velocity dispersion av- 
eraged over all galaxy types is cr = 160 km s^^, 
and the effective temperature is Ty = ^nipa^/k = 
1.8 X 10^ K. At plasma temperature Ty and our 
nominal heavy element abundance, Z — 10~^'^ 
times Solar, the Sutherland & Dopita^ (1993) cool- 
ing function is A 

characteristic cooling time at r2oo is 



10 ^'^ "^ erg cm'^ s ^, and the 



T200 = 



1.5n2ookTy 



An, 



e,200 



('^200 — "•e,20o) 



50 Gyr. (2) 



If the velocity dispersion a is independent of time 
the plasma density at r2oo in this expression scales 
as the square of Hubble's constant and the ratio of 



-'Here and in the rest of the paper we use the Suther- 
land & Dopita (1993) values for the mean mass per 
particle and the ratio of the electron to total particle 
number densities. We do not attempt to adjust the 
primeval helium abundance to the larger present estimate 
because the correction is small and the adjustment likely 
to be confusing. We use the corrected energy density in 
http://www.mso.anu.edu.au/ ralph/. 



cooling to Hubble times in a smooth corona varies 
with redshift as 

Ht2oo = 3.5[n„(l + zf + l- flm]-^/^- (3) 

Quasar heavy element absorption lines at 

projected distances < 100 kpc from galaxies 
(Churchill et al. 2000 and references therein) show 
that the coronae would have to have been enriched 
in heavy elements, but the abundance could be 
well below our nominal value. At primeval abun- 
dances the cooling time is twice equation (2). 

We conclude that the plasma cooling time near 
the virial radius of a present-day galaxy can be 
longer than the Hubble time, assuming what 
seems to be a reasonable heavy element abun- 
dance. This is what allows the development of a 
model for a massive corona that has not yet signif- 
icantly settled to condensed forms of baryons. The 
corona would have to have started accumulating 
at relatively low redshift. Equation (3) indicates 
that at redshift z = 3 the cooling time of a mas- 
sive halo at the smaller virial radius that obtained 
then is less than the expansion time then, even for 
a smooth halo. Thus we expect that the baryons 
that had collected at z = 3 largely collapsed to 
condensed forms within an expansion time. It is 
important that even at near primeval abundances 
the present cooling time in the outer part of a 
massive halo is not much greater than the Hub- 
ble time, meaning that only a modest degree of 
clumping of the plasma would have caused most 
of a massive corona to have dissipatively collapsed 
before the present epoch. It is reasonable there- 
fore that in some simulations 25 to 40% of the 
corona has ended up now as condensed matter 
(Dave et al. 2001). The puzzle is that the ob- 
served condensed mass fraction is much smaller, 
about 8% in the FP estimate. The resolution we 
are considering is that the corona was assembled 
in a smooth enough way that the cooling time 
for the bulk of the baryons in the galaxy once 
assembled was close to equation (2), assembly of 
the major pieces of a spiral galaxy having been 
largely completed at .2 > 1. 

3. The Corona Model 

We commence with a review of the FP conven- 
tions for the halo mass structure (Model 1) and 
then present an argument for an adjustment that 
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preserves the total mass but distributes it in a way 
that gives a better fit to the galaxy-mass cross cor- 
relation function (Model 2). Parameters appropri- 
ate for the massive halo of a spiral galaxy arc pre- 
sented in §3.2, the corona models are developed in 
§3.3, and the observational tests are discussed in 
§§3.3 and 3.4. 

3.1. The Mass Distribution 

In the halo model in FP, which is the basis for 
our corona Model 1, the total mass density p as a 
function of radius r is 



p{r) = 



27rGr2 ' 



r <ry = 



a 



50Ho 



(4) 



The circular velocity is constant at Vc = V^a at 
r <ry, the mass within the radius Vy is 



(5) 



and the mean density of mass within is 200 
times the Einstein-de Sitter value, that is, r« = 
r2oo- In the model to be discussed next we keep 
Tv at the value in equation (4) but adjust the mass 
distribution in a way that decreases r2oo- 

The FP estimate of the one-dimensional mass- 
weighted velocity dispersion averaged over galaxy 
types is 

(T = 160kms"^ (6) 
The effective number density of galaxies is 



Ug = 0.017/1^ Mpc" 



(7) 



where the Hubble parameter is Hq = 70 km 
Mpc^^, that is, h = 0.7. With these numbers the 
mass fraction within distance r„ from the center 
of a galaxy is 



f9 = 



MyUg 

^mPcrit 



0.6. 



(8) 



The baryon mass fraction 1 — fg beyond r„ might 
be assigned to the intergalactic medium, though a 
substantial part may be concentrated near galax- 
ies. 

In our corona Model 2 we adjust the halo mass 
distribution to produce a less abrupt truncation 
at r = r„ and, as will be discussed, a better ap- 
proximation to the galax;y-mass cross correlation 
function. One might consider smoothing the trun- 
cation of p(r) in equation (4) by adding mass at 




100 



radius, kpc 



1000 



Fig. 1. — Galaxy- mass cross correlation functions. 
The dashed line is the halo Model 1 in equa- 
tion (4), the solid line is the halo Model 2 in equa- 
tion (9) with r] = 1, and the dotted line is the 
SDSS weak lensing measurement. 



r > Ty with a steeper density gradient, but that 
would make the baryon mass in galaxies unreason- 
ably large, as one sees from the substantial fraction 
of mass within in equation (8) . We instead keep 
the effective galaxy mass fixed at My, thus hold- 
ing fixed the baryon fraction in equation (8), and 
spread out the mass distribution near ry. This is 
implemented in the model 



p{r) = 



2TrGr'^{l+r]r/ry) ' 
e'' - 1 

r<rx=ry , (9) 

V 

with p{r) = at r > rx- The total mass in this 

model agrees with My in equation (5) . The radius 
Ty here serves as an analog of the break radius a 
in the Hernquist (1990) form or in the Navarro, 
Frenk, & White (1996) model, but the form we 
are using has the steeper inner slope wanted to fit 
spiral rotation curves and weak lensing measure- 
ments. Note that the mass density near the center, 
set by (7, is the same as that in Model 1. The limit 
ry — > is equation (4). In our corona Model 2 we 
set T] = 1. 
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The model for the halo structure is tested by 
the galaxy- mass cross correlation function S^gp, 
which we approximate as (Seldner & Peebles 1977; 
Cooray & Sheth 2002) 



+- 



r2 f 
rv J 



2r^ J P{l+r]l/r^) 



(10) 



The step function, 9{r) = 1 at r < and 6{r) = 
at r > Tj., represents the extent of the halo. The 
first term on the right hand side is the contri- 
bution to the cross correlation function from the 
mass belonging to the galaxy. The second term is 
the convolution of the mass distributions in neigh- 
boring halos with the halo autocorrelation func- 
tion, which we approximate as ^ = (ro/r)^ with 
hvo = 5 Mpc and 7 = 1.8. The distance to a 
neighboring halo is d and I = -^/d^ — 2rd/i is 
the distance from the center of the neighbor to the 
point of observation of the mass density. This ex- 
pression assumes the mass fraction 1 — /g (eq. [8]) 
not in halos is homogeneously distributed, which 
seems reasonable at r ~ 100 kpc but at r > 1 Mpc 
may underestimate ^gp. 

The dashed line in Figure 1 is the galaxy-mass 
cross correlation function ^gp for the p oc 
halo model in equation (4), the solid line is S,gp 
for the halo model in equation (9), and the dot- 
ted line is the Sheldon et al. (2004) power law 
fit to their weak lensing measurement, ^gp = 
(5.4 Mpc//ir)^-^^. FP chose the p oc model to 
agree with the £,gp measurement as well as the cir- 
cular velocities observed at smaller radii. The fit 
to the former is indeed close, but ^gp is somewhat 
high at r ~ 200 kpc because halos overlap, and 
the abrupt cutoff in p{r) produces a pronounced 
shoulder in ^gp aX r ~ r^. The steeper halo den- 
sity run at r ~ in the 77 = 1 model removes 
the excess and reduces the shoulder. A still larger 
value for rj would make the shoulder even smaller 
but it would produce an undesirable dip in £^gp at 
r ~ 100 kpc. PAr > Tx both models fall below the 
measured ^gp. We attribute this to the mass frac- 
tion I — fg that we suspect is clustered on scales 

1 Mpc. 

We conclude that the models with 77 ^ and 
T] = 1 both give useful approximations to the mas- 
sive halos of L ~ galaxies, within the con- 
siderable uncertainties, but that the r] = 1 case 



produces the better fit to the galaxy-mass cross 
correlation function. 

3.2. The Massive Halo of a Spiral Galaxy 

We arc considering the idea that field galaxies 
have massive coronae with baryon mass fractions 
close to the cosmic mean. For late galaxy types 
it is appropriate to use a lower velocity dispersion 
than the average over types in equation (6). We 
adopt the FP value for spirals. 



CTiatc = 140 km s" 



(11) 



This defines the effective temperature, 

T„ = ^impalJk = 1.4 x lO^K, (12) 

and the late-type halo parameters (eqs. [4] and 
[5]), 

r„ = 280 kpc, My = 2.6 x (13) 

An estimate of the mean mass of the Milky Way 
and the Andromeda galaxy from the dynamics of 
relative motions of the nearby galaxies is 2.3 x 
IO^^Mq (Peebles 1995), reasonably close to M„ 
with our choice of criato- 

The conventional definition of the virialized 
part of a galaxy is that within the radius r^oo at 
which the enclosed density of mass is 200 times 
the Einstein-de Sitter value. In our late-type halo 
Model 1, with 77 — > in equation (9), r^oo = = 
280 kpc. In Model 2, with 77 = 1, the density 
of mass within is 40 times the Einstein-de Sit- 
ter value, or 140 times the present mean matter 
density. The lower density in the outer parts of 
this model reduces the conventional virial radius 
to r2oo = 240 kpc and it reduces the density pa- 
rameter of the mass within r2oo in spiral galax- 
ies from Q.y = 0.08 in Model I to Vty = 0.05 in 
Model 2. One should bear in mind, however, that 
the region of near dynamical equilibrium may ex- 
tend beyond r2oo- 

3.3. The Corona of a Spiral Galaxy 

In the FP inventory the condensed mass is 13% 
of the total baryon mass in a galaxy (from eq. [27] 
and entries 3.1a, 3.9 and 3.10 in Table 1 in FP), 
leaving the mass 



M,, 



O.SrMyQb/^m = 3.6 X 10"Mq (14) 
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for the massive corona of our characteristic spiral 
galaxy. 

Since the cooling time in the denser inner part 
of this corona is shorter than the Hubble time, and 
the Milky Way seems not to have been substan- 
tially disturbed by mergers since redshift z > 1 
(Gilmorc, Wyse & Norris 2001), we are assuming 
that the inner corona of a spiral typically has had 
time since the last major merger to have relaxed 
to a near steady rate of settling. This steady state 
means the particle number density n{r) and the 
inward streaming flow v{r) are such that the par- 
ticle flux entering the shell of radius r. 



N = A-KT nv, 



(15) 



is nearly independent of r outside the region where 
the condensed baryons are accumulating. The flux 
of entropy into the shell r is sN, where the entropy 
per particle is s = A; In T^/^/n, so the condition for 
a stationary energy distribution at radius r is 



kTN— In 

dr n 



= 4Trr^A{T)ne{n-ne). (16) 



Expressed in terms of the cooling time (as in 
eq. [2]) the particle number flux is 



N 



Gnnr 



dim 



d\nT^/^/n 



(17) 



Equations (15) and (16) are the steady-state ver- 
sions of the mass and energy equations used in 
analyses of cluster cooling flows (Fabian 1988; 
Sarazin 1988). 

Where the cooling time is greater than the 
free-fall time ^ r/a the plasma may be pressure- 
supported, and one might use the condition that N 
is constant with the equation of hydrostatic equi- 
librium, 

d ,^ GM(<r) 

—nekT= 5 fiUemp, (18) 

dr r"^ 

to compute the plasma temperature and density as 
functions of radius. We adopt a simpler approach, 
by noting that if the plasma temperature is inde- 
pendent of radius and the particle number density 
varies as n oc r~^/^ then equation (16) indicates 
that N is constant, as wanted under the condition 
that significant amounts of mass are dropping out 
of the corona only near the optical radius of the 
galaxy. 
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100 
radius, kpc 

Fig. 2. — Electron number densities for Model 1 
(the dashed line) and Model 2 (the solid line). 
Both take the plasma temperature to be T = 

4r^/3. 



In corona Model 1 we take the mass distribu- 
tion to be p oc (eq. [4]) and we take the plasma 
temperature to be Te = 4T„/3 so the plasma den- 
sity varies as r~^/^, making N constant. This is a 
special case of Cavaliere's (1973) beta model with 
/? = T^/Te = 3/4. The dashed line in Figure 2 
shows the electron number density in this model 
with the parameters in equations (11) to (14). At 
heavy element abundance 10~^-^ times Solar and 
plasma temperature Te = 4T„/3 = 1.8 x 10^ K 
the Sutherland-Dopita cooling fmiction is A = 
2Q-23.02 gj.g gj^3 g-1 j^^ primeval element abun- 
dance the cooling function is half this value. The 
mass flux dMcomi/dt = imipN in this model for 
the corona, computed from equation (16), is 



dM^ond/dt = UMq yr-i (lO"^-"^ Solar), 
dMcond/dt = 7Mq yr"^ (primeval). (19) 

This is the rate at which mass is entering the con- 
densed baryon components. It may be compared 
to the FP estimate of the present mean star for- 
mation rate per spiral galaxy (from the ratio of 
FP eqs. [30] and [52] with the spiral fraction 0.7), 



dM,t^,/dt = 1.7Mq y- 



(20) 
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Fig. 3. — The flux of mass out of the corona in 
Model 2 with heavy element abundance 10~^'^ 
times Solar. The upper curve is the ratio of the 
plasma cooling time to the Hubble time. The 
lower curve is the measure of the mass flux in equa- 
tion (16). Where the curve is dotted the cooling 
time is longer than the Hubble time and the ex- 
pression for the mass flux need not apply. 

This is less than the flux of mass out of the corona 
in Model 1 even at primeval element abundances. 
Since we know of no place where the condensed 
baryons might be accumulating at the rates in 
equation (19) the indication is that N in Model 1 
is too large by a modest but significant factor. 

We can remedy the problem by reducing rie- 
In Model 2 this is done by adopting the density 
profile in the modified halo model in equation (9) 
with rj = I while keeping the plasma temperature 
at T = 4T^/3. Since the mass distribution is not 
isothermal this is not a /3 model, but at r < r^, 
the mass density run is close enough to p ^ 
for our purpose. This is illustrated in Figure 3: at 
r < 60 kpc the cooling time is less than the Hubble 
time, so we want A'^ to be independent of radius, 
and we see that this is reasonably well satisfied at 
mass fiux 

dM^ond/dt ~ IMq yr-i (Model 2) (21) 

for heavy element abundance 10~^'^ times Solar, 
and half this value at primeval element abun- 



dances. At larger radii N shows a substantial vari- 
ation with radius but this need not matter because 
here the plasma cooling time is long and the steady 

flow condition need not apply. 

The mass flux out of the corona in Model 2 
(eq [21]) is satisfactorily close to the present star 
formation rate (eq. [20]). We continue the discus- 
sion of observational tests of the massive corona 
models in §3.4 and §4, after commenting on the 
nature of the plasma distribution in our models. 

In both corona models the plasma is warmer 
than the effective temperature of the dark matter. 
This moves the baryons toward the outer parts of 
the corona, making the baryon mass fraction in 
the central parts smaller than the cosmic value. 
This effect is more conspicuous in Model 2, where 
the halo is more extended, further lowering the 
electron density in the inner part of the halo by 
moving more baryons to the outskirts. The sup- 
pression of the baryon mass fraction in the inner 
parts of a massive halo is observed in clusters of 
galaxies, as shown by Allen et al. (2004) and ref- 
erences therein, and we assume the cause of the 
suppression may be the same for clusters and mas- 
sive galaxies: entropy is added to the plasma prior 
to and during assembly of the system. 

Figure 4 illustrates the situation in Model 2. 
The upper solid curve is the total mass contained 
within radius r (eq. [9] with 77 = 1) as a function of 
r. The lower solid curve is the sum of the corona 
and condensed baryon mass contained within r, 
and the dotted curves show the two separate com- 
ponents. We represent the condensed mass by the 
form 

Meo„d(< r) = (1 - e-'-/'-«) . (22) 

There is no intended physical significance to this 
functional form: we choose it to satisfy the con- 
dition that the condensed baryon mass density is 
close to the total at small radius and we fix the 
parameter Vg by the condition that the total con- 
densed mass, Mcond = '^'^iate'''g/^^ is 13 percent 
of the total baryon mass in the corona (as dis- 
cussed in §3.3). Under these conditions the ra- 
dius that contains half of the condensed mass is 
Te = Tg In 2 = 4 kpc;. This is in the range of effec- 
tive radii of L ~ L* spiral galaxies. Also, one sees 
in Figure 4 that in this representation the included 
condensed baryon mass at radius r = 10 kpc is 
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10 100 
radius, kpc 

Fig. 4. — The mass contained within a centered 
sphere in Model 2. The upper solid curve is the 
total included mass in baryons plus dark matter. 
The lower solid curve is the total baryon mass. 
The dotted curves show the two dominant baryon 
components, the plasma mass in the corona and 
the more concentrated baryon mass in condensed 
forms (eq. [22]). 

one half of the total, in agreement with what is 
known about the dark mass distribution in the 
Milky Way. Since the total condensed baryon 
mass agrees with the FP estimate from what is 
observed we conclude that this is a reasonable rep- 
resentation of the condensed mass distribution. 

At r = 20 kpc the included baryon mass is close 
to the total condensed mass, 6 x IO^^Mq, and the 
subdominant mass in the corona is causing the 
slow increase in included mass with increasing ra- 
dius. The included masses in condensed baryons 
and the corona plasma are equal at r = 180 kpc in 
Model 2. At this radius the total baryon mass frac- 
tion is 0.09, or 56 percent of the cosmic value. This 
suppression of the baryon fraction is comparable 
to what is observed in relaxed clusters of galax- 
ies at one quarter of the virial radius (Allen et al. 
2004). The significant difference is that in spirals 
only half of the baryons within r ^ 200 kpc are in 




5 10 15 20 

d sin |b|, kpc 

Fig. 5. — Dispersion measures along the lines of 
sight to pulsars in globular clusters. In order of 
increasing distance d sin \ b\ from the plane they are 
47 Tuc, M15, M13, M5, M30, NGC1851, M3, and 
M53. The lines assume a close to homogeneous 
corona with electron density rig = 0.001 cm~^ in 
the lower line and rig = 0.002 cm~^ in the upper 
line. 

the coronal plasma, while the other half is in con- 
densed forms. This difference might be expected 
from the shorter cooling times in galactic coronac. 
The included baryon mass fraction reaches the cos- 
mic mean at radius Vx (eq. [9]), at the right-hand 
edge of the figure. This would be replaced by a 
smooth transition to the cosmic value in a more 
realistic model for the corona. 

3.4. Observational Tests 

Parameters characteristic of our two models are 
summarized in Table 1. We have already discussed 
the entry for the rate of settling of the corona; the 
results for Model 1 (eq. [19]) motivate Model 2 
(eq. [21]). Here we consider other observational 
tests. 

If the Milky Way has a close to standard mas- 
sive corona then the dispersion measures of pul- 
sars outside the disk offer a useful probe of the 
plasma density in the inner corona. Manchester 
et al. (2005) list^ dispersion measures of pulsars 

^We use the catalog in 
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Table 1 
Massive Corona Parameters 





Model 1 


Model 2 


Electron density at r — 10 kpc, cm ^ 


0.010 


0.0026 


fie averaged from r — 10 to 50 kpc, cm~^ 


0.003 


0.0007 


Hqt at the outer edge of the halo 


2.6 


10.2 


o-latoT/r at r = 10 kpc 


4 


14 


Mass flux nN, Mq 


14 


1 


Compton-Thomson y parameter 


10-^" 


10-7.5 


Lo(e > 0.3 kcV, p > 10 kpc), erg 


^040.6 


1039.8 


Lo{e > 0.3 koV, 10 < p < 20 kpc), erg 


1040.0 


1038.8 


Soft X-ray background /gc, keV cm"^ st~^ keV~^ 


~ 1 


10 


Central photoionization rate, 


10-13.4 


10-14.6 



in the Magellanic clouds and in eight globular clus- 
ters that arc more than 2 kpc from the plane of 
the galaxy. Dispersion measures for the globular 
cluster pulsars are plotted in figure 5. The two 
lines assume a plane parallel plasma disk com- 
ponent with scale height <C 2 kpc and surface 
density 15 pc cm~^ together with a near homo- 
geneous corona. This value for the plasma surface 
density in the disk is not inconsistent with what 
is observed for pulsars closer to the central plane 
(Cordes & Lazio 2005, fig. 1). There is not com- 
pelling evidence for detection of a corona outside 
the disk component, but one sees from the fig- 
ure that the data would allow a corona with local 
number density in the range 

0.001 cm-3 <ne< 0.002 cm'^. (23) 

The electron density in Model 2 at r = 10 kpc, 
Ue = 0.0026 cm~"^, is reasonably consistent with 
this. In Model 1 rig is high, but we note in §4 
that this could be because clumping of the plasma 
substantially enhances the cooling rate and lowers 
Ue at r 10 kpc. 

Pulsars in the Magellanic Clouds have disper- 
sion measures ~ 100 cm~"^ pc, which we can take 
as a bound on the contribution from the corona. 
This corresponds to mean density 

ne(MC) ~ 0.002 cm-^ (24) 

along the line of sight to the Magellanic Clouds. 
Table 1 lists the mean electron density averaged 
in the radial direction from r = 10 kpc to r = 
50 kpc. The models are not inconsistent with 
equation (24). 

http: // www. atnf . csiro.au/ research / pulsar/psrcat 



Spitzer (1956) noted that the corona would be 
in pressure equilibrium with interstellar matter in 
the disk of the galaxy. Jenkins & Tripp (2001) 
find that the thermal pressure in the disk of the 
Milky Way within a few kpc radius averages nT = 
2000 cm""^ K, and they remark that this is not 
inconsistent with the weight of interstellar matter 
at the plane of the disk. In Model 2 the pressure at 
r = 10 kpc is nT = 5000 cm~^ K. We do not know 
how to judge whether the difference is significant. 

The integral of the electron density along a 
line of sight through the corona determines the 
effect of free electron scattering on the spectrum 
and anisotropy of the thermal 3 K background ra- 
diation (Rasmussen & Pedersen 2001 and refer- 
ences therein) . This is measured by the Compton- 
Thomson parameter y = kTo-T^e/ifneC^)- For 
the numbers listed in Table 1 we compute the 
surface density of electrons integrated radially 
from 10 kpc to the outer edge of the corona, ig- 
noring neighboring coronae. The empirical bound 
from the possible departure from a thermal spec- 
trum is 2/ < 1.5 X 10-5 (Fixsen et al.l996). The 
variation Sy of this parameter across the sky gives 
the contribution to the radiation anisotropy, 6y ~ 
ST/T <C 10"^, to avoid spoihng the consistency of 
the theory and observation of the anisotropy (Ben- 
nett et al. 2003). In both models y is well below 
these constraints. The same is true of the contri- 
bution of coronae to the mean electron pressure 
integrated across the Hubble length. 

The coronae we are considering are so large that 
the soft X-ray observation of a nearby galaxy con- 
strains the X-ray surface brightness rather than 
the total luminosity of the galaxy. In the usual 
representation in terms of an effective luminos- 
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ity, the integral of the surface brightness over 
the observed part of the galaxy is a recieved en- 
ergy flux density / and the effective luminosity is 
Lq = 47rD^/ if the distance to the galaxy is D. 
We compute this luminosity derived from a mea- 
surement of the energy flux density received from 
a cylinder of radius p centered on the galaxy as 



Lq = An 



f 

Jo 



r^drF{r)j{r), 



(25) 



where F{r) is the fraction of a sphere of radius r 
centered on the galaxy that falls within the cylin- 
der of radius p. Thus F = 1 at r < p and 



F{r) = 1 - a/1 - p2/r2 



(26) 



at larger radii. We approximate the plasma X- 
ray luminosity density at energy greater than e as 
(Spitzer 1978) 

j = 1.7X 10-27n2TV2e-^/'=^ erg cm'^ s-\ (27) 

in cgs units with the plasma temperature in °K. 

The X-ray luminosities in Table 1 are computed 
at energy e > 0.3 kcV, to match the measurements 
by Strickland et al. (2004a), and for the two mod- 
els of the corona. The first of the entries refers to 
the energy flux density integrated from projected 
distance p= 10 kpc to p = at the outer edge of 
the corona. The second entry refers to the energy 
flux density from the range of projected distances 
10 < p < 20 kpc, which is closer to available ob- 
servations. There are several reasons to exclude 
the part at p < 10 kpc: neutral hydrogen strongly 
absorbs soft X-rays, the corona plasma likely is 
strongly perturbed by the disk, and X-ray sources 
in the disk arc not readily separated from what 
might be produced by the corona. 

Strickland et al. (2004a) flnd that the halo 
regions of the nearby galaxies NGC 891 and 
NGC 253, both of which have circular velocities 
similar similar to the Milky Way, have luminosi- 
ties ~ 10"^^ erg s"'^ at 0.3 to 2 keV. This is com- 
parable to the luminosity at 10 < p < 20 kpc in 
Model 2 and an order of magnitude below Model 1. 
The significance is unclear, however, because what 
is observed may have been produced by matter 
blown out of the disks by present or past star- 
bursts, rather than the effect of a corona (e.g. 
Pietsch et al. 2000; Strickland et al. 2004b), and 



one must bear in mind that the luminosity com- 
puted for Model 1 may be an overestimate for the 
reason discussed in §4. 

The total corona X-ray luminosity is con- 
strained by the possible contribution to the soft 
X-ray background, which Wu, Fabian, & Nulsen 
(2001) put at 



/ < 4 keV cm-2 s"^ sr'^ keV" 



(28) 



at photon energy e ~ 0.25 keV, after subtraction 
of detected sources. The contribution to / by coro- 
nae depends on their history, which is even more 
uncertain than their present state, so we use a 
crude approximation. 



(29) 



The characteristic galaxy number density is Ug 
(eq. [7]), the spiral fraction is fs = 0.7, we approx- 
imate the time integral as J L dt/ (1 + z) ~ O.bLt, 
with t = 10 Gyr and L constant at the values in 
Table 1 for p > 10 kpc, and the factor e = 0.3 keV 
in the denominator approximatc!s the width of the 
spectral energy distribution. This estimate, listed 
in the second line from the bottom of Table 1, 
ignores the shift of the photon energies out of 
the band of detection, so if L were close to con- 
stant it would be a considerable overestimate. Wu, 
Fabian, & Nulsen (2001) point out that the corona 
X-ray luminosity could have been substantially 
larger in the past, meaning equation (29) could 
be an underestimate, but that depends on how 
the galaxies and their coronae formed and evolved. 
We conclude that our picture for the present-day 
state of massive coronae in spiral galaxies does not 
conflict with the soft X-ray background measure- 
ment, within the considerable uncertainties. 

The X-ray emission from coronae contributes to 
the local ionizing radiation, as Maloney & Bland- 
Hawthorn (1999) note in their discussion of plasma 
distributed through the Local Group. An observer 
situated near the center of a corona would find 
that the radiation energy density produced by the 
corona is roughly isotropic at 



1 



j{r)dr, 



(30) 



10 kpc 



where j{r) is the luminosity density (eq. [27]). As 
before, we assign the lower limit r = 10 kpc be- 
cause the situation at smaller radii is likely to 
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be dominated by the disk. With photoionization 
cross section a oc v~'^ ^ and a{yo) = lO"^''-^ cm^ at 
the hydrogen photoionization threshold, the pho- 
toionization rate near the center of the corona is 

r ~ (ToCu/(3A;T), (31) 

The numerical results arc listed in Tabic 1. The 
Penton et al. ( 2000) estimate of the rate of ion- 
ization by the intergalactic background radiation, 
r = Vd~^^-^ s"-'^, is comparable to what might be 
expected from the corona in Model 1. In Model 2 
it appears that massive coronae would not be sig- 
nificant local sources of ionizing radiation. 

4. Discussion 

We have argued for the viability of the idea 
that about half of the baryons are in gravitation- 
ally bound coronae in the dark matter halos of 
galaxies. In our models the mass distribution is 
designed to fit the galaxy-mass cross correlation 
function (Fig. 1), under the assumption of near 
monolithic massive halos, and the coronae are as- 
sumed to have relaxed to a steady cooling flow. 
Under these conditions the corona electron den- 
sity run is not likely to be much outside the two 
curves in Figure 2: a lower density than in Model 2 
would bring the density in the outer parts of the 
corona closer to the global mean than seems rea- 
sonable for matter that has relaxed to dynamical 
equilibrium, and a density greater than or compa- 
rable to Model 1 would cause substantial problems 
with the observations. 

A better massive halo model would replace the 
sharp truncation of the mass and baryon densi- 
ties at r = with a physically realistic tran- 
sition to the intergalactic medium, but what is 
relevant for the inner structure of a corona in a 
near steady state is the situation where the plasma 
cooling time is comparable to the Hubble time, at 
r ~ 60 kpc: the plasma with longer cooling times 
serves as a reservoir supplying the flow to smaller 
radii. Conditions at this transition to near steady 
flow likely depend on departures from spherical 
symmetry, as may be caused by tidal interactions 
with neighboring massive halos or ongoing accre- 
tion. One might also have to consider rotational 
support in the inner parts of the corona, and possi- 
bly magnetic field pressure. But such refinements 
might await clearer evidence that massive coronae 
really exist. 



Model 2 appears to be observationally accept- 
able. It should be noted, however, that this model 
depends on the assumption of a smoothly dis- 
tributed plasma, which may not be accurate, be- 
cause the plasma is thermally unstable to the 
growth of density fluctuations that decrease the 
cooling times and so reduce the electron density 
run from what is shown in Figure 1.'^ That may 
make the X-ray luminosity of the corona signifi- 
cantly smaller than our estimate. The thermal in- 
stability may be expected to have a lesser effect on 
the flux of mass from the corona into condensed 
forms of baryons, because this is determined by 
conditions where the cooling time and the ther- 
mal instability growth time both are comparable 
to the Hubble time. 

Model 1 appears to be marginally unaccept- 
able: the most serious problem is the large rate 
of production of condensed baryons (cq. [19]). It 
is worth noting, however, that if galactic winds 
removed about half of the baryons in the corona, 
leaving the rest with the radial distribution in 
Model 1, it would reduce the mass fltix by a factor 
of four, leaving a still substantial baryon mass in 
coronae and an acceptable and interesting present 
rate of production of condensed forms of baryons. 

The near coincidence between the flux of mass 
out of the corona (eq. [21]) and the present star 
formation rate (eq. [20]) in Model 2 — or in 
Model 1 with half the plasma mass we have been 
assuming — agrees with Larson's (1972) con- 
jecture that gas is flowing into the interstellar 
medium at the rate the interstellar medium is be- 
ing converted to stars. Inflow to the disk at about 
this rate figures in current discussions of the his- 
tory of the thin disk of the Milky Way (as reviewed 
by Pagcl 1997; Matteucci 2001). The high veloc- 
ity clouds often are mentioned as a source for the 
inflow, but the absence of analogs around other 
galaxies makes this unlikely (Pisano et al. 2004). 
Mass settling out of the corona seems to be a vi- 
able alternative. In this picture star formation in 
an isolated spiral galaxy will continue for a con- 
siderable time to come. 

We have not analyzed compatibility of the mas- 

^In the range of radii 10 < r < 60 kpc the corona coohng 
time is less than the Hubble time but larger than the free 
fall time. Here we may apply the quasi-static approxima- 
tion in eq. [17], where the density scales in proportion to 
the cooling time. This scaling ends at free fall. 
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sive corona picture with the much cooler clouds of 
atomic hydrogen and low ionization ions detected 
near galaxies by quasar absorption line observa- 
tions. Perhaps the plasma thermal instability pro- 
duces the cloudy multiphase structure indicated 
by quasar absorption line studies (Churchill et al. 
2000 and references therein). Chen et al. (2001) 
find that the HI surface density at projected dis- 
tance p = 100 kpc from an L, galaxy at modest 
redshift is typically ^ 10^"^ cm~^. Since the sur- 
face density of free electrons at p = 100 kpc is 
10^°'^ cm~^ in Model 2 there would be enough 
mass to produce the observed Lya absorption. 
The minimum size of a region dense and cool 
enough to allow significant densities of HI and low 
ionization ions likely is set by thermal conduction 
by the motions of hydrogen atoms between forma- 
tion and ionization. But we have left for future 
work the subtle analysis of the surface densities of 
these cool clouds. 

Another subject for more careful discussion is 
the history of the flux of mass out of the coronae 
into condensed forms of baryons. The star forma- 
tion rate has to have been larger in the past, as has 
been widely discussed: an example of the evidence 
is the difference between the present star forma- 
tion rate in equation (20) and the mean rate of 
production of the observed present mass in stars 
and interstellar matter averaged over the age of 
the universe, 

(dM,,.,/dt) = ^h^^^ = y-\ (32) 

UgtQ 

where Ocoid = 0.0035 is the FP estimate of the 
condensed baryon mass and Iq = 14 Gyr. Within 
the massive corona picture it is easy to imagine 
why the flux of mass out of the corona and into 
stars might be decreasing. We have estimated that 
at redshift 2 = 3 the cooling time in the virialized 
mass concentration around a galaxy is less than 
the expansion time (eq. [3]). That is, a corona as- 
sembled then would collapse to condensed matter 
in the disk or bulge. At lower redshift hierarchical 
growth of a galaxy would produce density fluc- 
tuations that would tend to increase the cooling 
rate and hence increase the flux of mass out of a 
corona present then. Settling of the denser clumps 
of plasma would tend to leave behind a smoother 
longer-lived corona. And at low redshifts super- 
novae may add entropy to the corona, further sup- 



pressing the rate of dissipative settling. 

We showed in §2 that the massive corona pic- 
ture requires a close to smooth plasma distribu- 
tion, not seriously disturbed by merging and ac- 
cretion at low redshifts. There is empirical sup- 
port for this condition, from the evidence that 
the Milky Way has not suffered a "major inva- 
sion" since redshift unity (Gilmore, Wyse & Nor- 
ris 2001), and from the absence of the HI clouds 
around galaxies (Pisano et al. 2004) that would 
seem to be a natural consequence of merging and 
accretion at low redshifts (Blitz et al. 1999). 

The condition that we can neglect the effects 
of merging and accretion is not suggested by nu- 
merical simulations of the ACDM cosmology (as 
illustrated in Fig. 2 in Gao et al. 2004). It may 
be important therefore that a possible remedy is 
under discussion: a long-range scalar interaction 
in the dark sector would hasten the assembly of 
dark matter halos and thus suppress the rate of 
major halo mergers at 2; < 1 (Nusser et al. 2005 
and references therein). Since baryons cannot sig- 
niflcantly couple to the scalar force in this pic- 
ture they would tend to accumulate in halos more 
slowly than the dark matter. This suggests yet 
another picture for the disposition of the baryons: 
maybe a substantial fraction have not fallen into 
massive halos. Or perhaps the galaxies did collect 
the full baryon mass fraction, but in a gentle way 
that allowed the formation of long-lived massive 
coronae. 

This paper is motivated by a challenging issue 
— where are the baryons that are expected to have 
accompanied the dark matter in the formation of 
a galaxy? — and by a natural idea - that most 
of these baryons are in pressure-supported plasma, 
the coronae, with radii comparable to that of the 
dark matter halos of L ^ galaxies. We con- 
clude that the idea is viable but needs work. On 
the observational side, we will be following with 
particular interest the development of measure- 
ments of the soft X-ray surface brightnesses of iso- 
lated spiral galaxies at projected distances in the 
range of 100 to 300 kpc from the galaxy, which 
we expect would sample the bulk of a corona and, 
at the natural corona temperature in the range of 
one to a few times 10® K, would be detectable if 
present. On the theoretical side, pressing issues 
are the compatibility of the corona picture with 
the cooler clouds identified by quasar absorption 



12 



line studies and compatibility with the theory of 
structure formation within some acceptable ver- 
sion of the standard cosmology. 
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